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Abstract 



> , 

. AO 0235+164 is a very compact, flat spectrum radio source identified as a BL Lac object at a redshift of 

z = 0.94. It is one of the most violently variable extragalactic objects at both optical and radio wavelengths. 
The radio structure of the source revealed by various ground-based VLBI observations is dominated by a 
' nearly unresolved compact component at almost all available frequencies. 
. Dual-frequency space VLBI observations of AO 0235-fl64 were made with the VSOP mission in January- 

' February 1999. The array of the Japanese HALCA satellite and co-observing ground radio telescopes in 
r~| . Austraha, Japan, China and South Africa allowed us to study AO 0235-1-164 with an unprecedented angular 
Qlh' resolution at frequencies of 1.6 and 5 GHz. We report on the sub-milliarcsecond structural properties of 
the source. The 5-GHz observations led to an estimate of Tq > 5.8 x 10^'^ K for the rest-frame brightness 
^ ' temperature of the core, which is the highest value measured with VSOP to date. 

C/3 i Key words: Galaxies: active — Radio continuum: galaxies — Radio sources: variable — BL Lacertae 

' objects: individual (AO 0235+164) — Techniques: interfcrometric 



5_] ■ 1. Introduction 1992), optical (Rabbette et al. 1996) and X-ray {Ma.de- 

[ jski et al. 1996) regimes of the electromagnetic spec- 

The Arecibo Occultation source AO 0235+164 is a trum. It is detected in gamma-rays by the EGRET on 
well known object extensively studied at nearly all wave- the Compton Gamma- Ray Observatory (von Montigny et 
lengths. It was identified as a BL Lac object by Spinrad al. 1995), with indication of variabihty during the 1-year 
and Smith (1975) based on its almost featureless opti- observing period. The presence of a group of foreground 
cal spectrum. Its emission line redshift of z = 0.940 galaxies led Stickel, Fried & Kiihr (1988) to suggest that 
was first measured by Cohen et al. (1987). Two other the dramatic flux density variations, correlated at opti- 
absorption redshift systems were detected in the same cal and radio wavelengths, could be explained in terms 
direction at z = 0.851 and z = 0.524. The latter was de- of gravitational microlensing. Abraham et al. (1993) 
tected also in emission. There are a number of foreground performed deep imaging of AO 0235+164. Based on the 
galaxies within a few arsceconds of AO 0235+164 at the observation and modcUing its Mgii absorber, they con- 
latter redshift, indicating the presence of a group or clus- eluded that strong microlensing of AO 0235+164 by in- 
ter of star-forming galaxies. Extensive optical studies of dividual stars in the intervening galaxy is unlikely. Saust 
the object AO 0235+164 and the surrounding field, both (1992) came to a similar conclusion using emission line 
imaging and spectroscopy, are presented in Nilsson et al. variability measurements. It appears that variability is 
(1996), Burbidge at al. (1996) and references therein. intrinsic to the source and may be best explained by the 

The source is highly variable in the radio (O'DeU et al. relativistic beaming model. 

1988; Romero et al. 1997), near- infrared (Takalo et al. Flux density variations across a factor-of-50 range 
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Table 1. Imaging VLBI observations of AO 0235+164 from the literature. 
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flux density is the sum of the fitted model component flux densities cited in the VLBI reference (where available) 
^ flux density of the core 

Gaussian model size of the core (FWHM) 

lower limit for the core brightness temperature estimated from fltted model flux densities and angular sizes 
" EVN=European VLBI Network; SVLBI^space VLBI; VLBA=Very Long Baseline Array 

^ AG95=Ahschuler et al. 1995; BE81=Baath et al. 1981; CB96=Chu et al. 1996; CU83=Cohen et al. 1983; CZ99=Chen, 
Zhang & Sjouwerman 1999; FFOO=Fomalont et al. 2000; GC89=Gabuzda et al. 1989; GC92=Gabuzda et al. 1992; 
GC96=Gabuzda & Cawthorne 1996; JB84=Jones et al. 1984; KV98=Kell ermann et al. 1998; SW97=Shen et al. 1997; 
USNO=US Naval Observatory ( ^ttp:/ /maia. usno.navy. mil/ rorf/ rrfid. htm^ ; VOO=T. Venturi et al. 2000, in preparation; 
XW98=Xu, Wehrle & Marscher 1998 

^ peak brightness and restoring beam size in the VLBI image were used to substitute S and size, to obtain a rough 
estimate of Tb, assuming that the core is unresolved 



in radio frequencies during a 5-year monitoring interval 
were found to be correlated, suggesting that they are 
intrinsic to the object. The radio variabihty can qual- 
itatively be understood in terms of adiabatically evolv- 



ing structures in a relativistic jet (O'Dell et al. 1988). 
However, Romero et al. (1997) favor refractive scintilla- 
tion as the simplest explanation of intraday radio vari- 
ability observed in AO 0235+164. Most recently, Kraus 
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et al. (1999) analysed contemporaneous radio and opti- 
cal monitoring data taken in October 1992 over a pe- 
riod of about a month. They discuss several models 
to explain the observations, including relativistic shock- 
fronts, processing beams, free-free absorption by the fore- 
ground medium, interstellar scattering and gravitational 
microlensing. Since apparently none of the models inves- 
tigated accounts for the observed variations alone, Kraus 
et al. suggest that the variability is caused by a super- 
position of intrinsic and propagation effects. They note 
that all explanations imply a very small intrinsic source 
size and require a Doppler factor substantially larger than 
usual (up to 100) to be consistent with the inverse Comp- 
ton brightness temperature limit of ~ 10^^ K (Keller- 
mann & Pauhny-Toth, 1969). 

A large number of VLBI images of AO 0235+164 have 
been produced over a wide range of observing wave- 
lengths (from millimeters to meters). We have collected 
a long, albeit certainly incomplete, list of references in 
Table 1. A large fraction of the VLBI images do not 
show any extended structure apart from the very com- 
pact core, regardless of the different angular resolution 
and the actual total flux density at the observing time 
(e.g. Gabuzda et al. 1989, 1992, 1996; Xu, Wehrle & 
Marscher 1998). On the other hand, several authors re- 
port faint extensions at a variety of position angles (PA) , 
but mainly between the north and cast (e.g. Jones et al. 
1984; Chu et al. 1996; Shcn et al. 1997). Although the 
components foimd by Chu et al. (1996) appear to move 
too fast to be identifiable between the subsequent epochs, 
these authors, based on the assumption that the compo- 
nents they observe are the results of outbursts at earlier 
epochs, come to an estimate of the Doppler boosting of 
the order 10^. 

Several non-imaging VLBI surveys, e.g. Preston et 
al. (1985) at 2.3 GHz, Morabito et al. (1986) at 2.3 
and 8.4 GHz and Moellcnbrock et al. (1996) at 22 GHz, 
also fomid AO 0235+164 to be extremely compact. The 
source is often used as a fringe-finder and calibrator for 
VLBI experiments. 

2. VSOP space VLBI observations, calibration 
and data reduction 

AO 0235+164 was observed in the VLBI Space Obser- 
vatory Programme (VSOP, Hirabayashi et al. 1998) at 
1.6 GHz on 31 January 1999 with the HALCA satellite 
via the tracking stations at Tidbinbilla (Australia), Gold- 
stone (USA), Usuda (Japan) and Robledo (Spain). The 
participating ground telescopes were the Australia Tele- 
scope Compact Array (ATCA), Sheshan (China), Har- 
tebeesthoek (South Africa) and Noto (Italy). The 5- 
GHz VSOP observation took place on 1 February 1999. 
HALCA was tracked from Tidbinbilla. The three co- 
observing ground radio telescopes were ATCA, Sheshan 



and Usuda. The data in both experiments were recorded 

in left circular polarization in S2 format (Wietfeldt et 
al. 1996) with 32 MHz bandwidth. Correlation was per- 
formed at the Dominion Radio Astrophysical Observa- 
tory in Penticton (Canada) (Carlson et al. 1999). 

Data calibration and fringe-fitting with 5-min solution 
interval were done with the NRAO AIPS package (version 
150CT98; e.g. Cotton 1995, Diamond 1995). Time aver- 
aging, editing, self-calibration and imaging were done us- 
ing the DIFMAP program (version 2.3c; Shepherd, Pear- 
son & Taylor 1994). 

Strong fringes were found in the 5- GHz experiment 
for all available antennas and time ranges. Measured 
system temperatures were available at Sheshan and 
Usuda for initial amplitude calibration. Nominal val- 
ues were used for HALCA and ATCA. The calibra- 
tor source 0420—014 was used to adjust the antenna 
gains of the ground telescopes. The total flux den- 
sity of 0420—014 is being monitored at the University 
of Michigan Radio Astronomy Observatory (UMRAO, 
http://wwwMstroAsaMmich.edv/obs/radiotel/umrao.html). 
The source was observed in the 5-GHz VSOP-VLBA Pre- 
launch Survey (VLBApls, June 1996, Fomalont et al. 
2000). Based on the similar 5-GHz total flux densities 
at the epochs of the VLBApls and our VSOP observa- 
tion, and assuming that the source structure found in the 
VLBApls has not changed considerably since 1996, the 
ground telescope gains had to be corrected by factors of 
1.3, 1.5 and 3.3 for ATCA, Sheshan and Usuda, respec- 
tively. The unusually large correction factor for Usuda 
is explained by the fact that the quantizer threshold was 
set incorrectly for this telescope, as noted at the correla- 
tor. Initial HALCA amplitude calibration could not be 
checked but it is known to be stable (Moellenbrock et al. 
2000). 

The (w, v) coverage for the 5-GHz experiment is shown 
in Fig. 1. Fig. 2 shows the visibility amplitudes and 
phases as a function of projected baseline length up to 
420 million wavelengths (MA), along with the curves rep- 
resenting the final clean-component model. The uni- 
formly weighted image is shown in Fig. 3. 

In the case of the 1.6- GHz experiment, the initial am- 
plitude calibration was based on the measured system 
temperatures and gain curves for Hartebeesthoek, Noto 
and Sheshan. Nominal system temperatures were used 
for ATCA and HALCA. Fringe-fitting resulted in strong 
fringes for all available baselines except those to She- 
shan, where no fringes were found. This, and the fact 
that the source was observed during a very limited com- 
mon time range at the remaining antennas, made the 
traditional hybrid mapping of AO 0235+164 impossible 
at this observing frequency. Model fitting to the visibil- 
ity data in DIFMAP resulted in an elliptical Gaussian 
component with flux density of 1.44 Jy, angular size of 
0.85 mas x 0.60 mas (FWHM) and major axis position 




Fig. 1.. The {u,v) coverage of the 5-GHz space VLBI experiment. 
Participating ground antennas are ATCA, Sheshan and Usuda. 
Projected baselines to HALCA extend up to 420 MA. 
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Fig. 2.. Self-calibrated visibility amplitude (Jy) and phase (de- 
grees) vs. projected baseline length for AO 0235-1-164 at 5 GHz. 
Solid lines represent the clean component model used for restor- 
ing the image in Fig. 3. 



o - 

o ; 

o - 



Fig, 



40 60 80 100 120 140 

UV rodius (10^ A) 

4.. Calibrated visibility amplitude (Jy) and phase (degrees) 
vs. projected baseline length for AO 0235+164 at 1.6 GHz. 
Reliable amplitude self-calibration could not be made due to 
the lack of four useful antennas in the array for most of the 
observing time. 



angle of —14°. Fig. 4 shows the visibility amplitudes 
and phases vs. projected baseline length. The data were 
edited and averaged over 30 s intervals in DIFMAP. The 
ground telescope antenna gains were again adjusted us- 
ing the total flux density data taken simultaneously and 
assuming that the source is unresolved on the shortest 
baselines. The gain correction factors were 1.4, 1.4 and 
1.8 for Noto, Hartebeesthoek and ATCA, respectively. 



3. Discussion 
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3.1. Sub-milliarcsecond structure 

At 5 GHz, the structure of AO 0235+164 on sub- 
milliarcsecond scale is dominated by a compact compo- 
nent, with a weak extension to approximately the north- 
west (Fig. 3). Its position angle is not quite the same 
as that found in many VLBI images pubUshed earlier 
(see Table 1), including the 22-GHz ground VLBI im- 
age of Jones ct al. (1984) which has ahnost exactly the 
same angular resolution. Note that earlier ground-only 
VLBI imaging observations led to controversial results 
on the existence and position angle of the extended fea- 
tures, which most likely indicates rapid structural varia- 
tions. However, the existence of such an extension must 
be treated with some caution because of the relatively 
sparse {u,v) coverage in our space VLBI observation. 
The visibility data of AO 0235+164 are consistent with 
a central point-source and an elliptical Gaussian compo- 
nent, according to model fitting in DIFMAP. The model 
parameters are listed in Table 2. Indeed, the source is 
remarkably compact: the correlated flux density ratio 
between the longest and shortest projected baselines is 
3:5 (Fig. 2). 

3.2. Source spectrum 

At the time of our observations, the total flux density 
of AO 0235+164 was declining according to the data from 
the UMRAO long-term monitoring program at three dif- 
ferent radio frequencies (Aller et al. 1999). The total flux 
density data at 4.8 GHz are shown in Fig. 5, indicating 
the epoch of our observations. The quasi-single epoch 
source spectrum (Fig. 6) is constructed from measure- 
ments at UMRAO and the Arecibo Observatory on 30 
January 1999. The flux density over a range spanning an 
order of magnitude in frequency (from 1.175 to 14.5 GHz) 
is roughly constant. It can be well described with a 
power-law spectral index, a = —0.02 ± 0.03 {S ~ v"). 
Although there are no measurements available at lower 
frequencies, a spectral turnover at ~ 2 GHz may be sus- 
pected. 

It should be noted that the relatively large ground ra- 
dio telescope gain correction factors used for amplitude 
calibration of our space VLBI data lead to correlated flux 
densities of ~ 1.5 Jy on the shortest baselines. This is in 
agreement with the contemporaneously measured total 
flux densities and the well-known compact nature of AO 
0235+164 which makes this source a widely used VLBI 
fringe-finder and calibrator. 

3.3. Brightness temperature 

A lower limit for the rest-frame brightness tempera- 
ture of the source can be estimated as 



4.8 GHz 




time (year) 



Fig. 5.. Total flux density of AO 0235+164 at 4.8 GHz measured 
at the University of Michigan Radio Observatory (UMRAO) 
26-m antenna since 1979. The vertical line indicates the epoch 
of our observations. 



0235+164 on 30 Jan 1999 



■^1 2 5 10 

frequency (GHz) 

Fig. 6.. Total flux density of AO 0235-1-164 at nine different fre- 
quencies measured nearly contemporaneously at the Arecibo 
Observatory (open circles) and at the UMRAO (triangles). Val- 
ues at lower frequencies are taken on 30 January 1999. Flux 

densities at 8 and 14.5 GHz arc linearly interpolated using the 
closest in time UMRAO data. 
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Table 2. Fitted model parameters of the source components at 5 GHz. Angular size upper limits and uncertainties are 
calculated from the image parameters, following Fomalont et al. (2000). 



Component 


S 


r 


e 


a 


6 






(mJy) 


(mas) 


n 


(mas) 


(mas) 


n 


Core 


915 


0.04 


127 


< 0.05 


< 0.03 




A 


319 


0.21 


-57 


0.60 ± 0.02 


0.20 ± 0.02 


37 



S: flux density, r: separation from the fleld center, 6: position angle, a: major axis, b: minor axis, position angle of 
the major axis; position angles are measured from north through east. 



where S is the component flux density in Jy, 0maj and 
(^min are the major and minor axes, respectively, in mas 
(assuming optically thick Gaussian brightness distribu- 
tion), u is the observing frequency in GHz and z is the 
redshift. For the core of AO 0235+164, at the redshift of 
z = 0.94, we obtain Tb > 5.8 x 10"^^ K using the source 
model from the 5-GHz experiment (Table 2). This is 
a direct evidence for a brightness temperature exceed- 
ing the inverse Compton limit and can be interpreted, in 
particular, as a result of relativistic motion in the source. 

Earlier brightness temperature estimates for AO 
0235+164 were based on e.g. 22-GHz VLBI survey data 
(Tb > 2 X 10^2 K, Moellenbrock et al. 1996), 5-GHz 
ground-based VLBI imaging (Tb = 5.3 x lO^^ K, Shen et 
al. 1997) and variability (Tb = 7 x 10^^ K for A = 20 cm, 
Kraus et al. 1999). We have made an attempt to com- 
pile a list of brightness temperature lower limit estimates 
using the VLBI imaging experiments cited in Table 1 , ac- 
cording to Eq. 1. Calculated Tb values are based on fitted 
Gaussian core model component flux densities and angu- 
lar sizes given in the references where available. The 
brightness temperature was estimated using image pa- 
rameters (peak brightness and restoring beam size, as- 
suming that the source is not resolved) where only images 
were available. The values as a function of time are plot- 
ted in Fig. 7 using the same time scale as for the total flux 
density monitoring data (Fig. 5). The estimates of Tb arc 
lower limits, the time sampling is quite sparse, and the 
values were calculated from VLBI experiments made at 
different frequencies, and thus with different sensitivities, 
to determine the brightness temperature. Nevertheless, 
there is an indication of a correlation between the long 
term behaviour of Tb and measured total flux densities 
(Fig. 8). The correlation cannot arise from the depen- 
dence of Tb on S alone, since source size is also highly 
variable (see Table 1) and Tb varies in a range spanning 
four orders of magnitude (Fig. 7). This possible cor- 
relation is qualitatively consistent with the picture that 
each flux density outburst marks the birth of a new, rela- 
tivistic component which later fades away and the angle 
between the direction of motion and the line of sight in- 
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1980 1985 1990 1995 2000 

time (year) 

Fig. 7.. Estimated lower limits of brightness temperature of 
AO 0235-1-164 vs. time. Open circles and filled boxes indi- 
cate estimates based on fitted VLBI model components and 
image parameters, respectively (see footnotes to Table 1 for 
explanation). An estimate of Moellenbrock et al. (1996) using 
non-imaging VLBI survey data is shown with a cross. Times 
of local maxima and minima in the 4.8-GHz total flux density 
curve (Fig. 5) are indicated with large tick marks at the top 
and bottom of the plot, respectively. 

creases (Chu et al. 1996). As the jet direction in BL 
Lac objects is thought to be oriented very close to the 
line of sight, this change causes significant decrease in 
the Doppler boosting factor. 

4. Conclusions 

Our 5-GHz VSOP space VLBI experiment shows the 
BL Lac object AO 0235+164 to be very compact on sub- 
milliarcsecond angular scales. The rest-frame brightness 
temperatures estimated from the source angular size and 
observed flux density clearly exceed 10^^ K which can be 
explained by strong Doppler boosting. The brightness 
temperature Tb > 5.8 x 10^^ K obtained from the 5-GHz 
data is the highest value found with VSOP observations 
to date (cf. Lovell et al. 2000). This implies a Doppler 
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1980 1985 1990 1995 2000 

time (year) 

Fig. 8.. Standardised values of the 4.8-GHz total flux density (So, 
dots; Fig. 5) and of the logarithms of brightness temperature 
lower limits (Tb, circles; Fig. 7) arc plotted together as a func- 
tion of time, to illustrate the apparent correlation between the 
two data sets. The correlation coefficient calculated using the 
closest-in-time So and Tb values is 0.7. 

factor of ^ 100, in good agreement with the results of 
several recent studies (e.g. Kraus et al. 1999, Fujisawa 

ct al. 1999), although the source total flux density was 
relatively low at the time of our observations. 

There is an indication of an extended component to 
N-NW, at the position angle about 90° away from what 
is sometimes scon in ground VLBI images. Due to the 
possibly very small jet angle with respect to the line of 
sight (e.g. (f) < 0?35 is estimated by Fujisawa ct al. 1999), 
small deviations in the intrinsic direction of the jet ejec- 
tion could cause large changes in the observed jet position 
angle. 

Further high resolution VLBI monitoring observations, 
either with space VLBI and/or ground-based VLBI at 
higher frequencies {u > 15 GHz), would be needed to 
identify possible components ejected after major out- 
bursts occur ing every 3-5 years. 
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